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ABSTRACT 

We report the detection of 24 /im variations from the planet-hosting v Andromedae system consistent 
with the orbital periodicity of the system's innermost planet, v And b. We find a peak-to- valley phase 
curve amplitude of 0.00130 times the mean system flux. Using a simple model with two hemispheres of 
constant surface brightness and assuming a planetary radius of 1.3 Rj gives a planetary temperature 
contrast of > 900 K and an orbital inclination of > 28°. We further report the largest phase offset 
yet observed for an extrasolar planet: the flux maximum occurs ~ 80° before phase 0.5. Such a large 
phase offset is difficult to reconcile with most current atmospheric circulation models. We improve on 
earlier observations of this system in several important ways: (1) observations of a flux calibrator star 
demonstrate the MIPS detector is stable to 10 -4 on long timescales, (2) we note that the background 
light varies systematically due to spacecraft operations, precluding use of this background as a flux 
calibrator (stellar flux measured above the background is not similarly affected), and (3) we calibrate 
for flux variability correlated with motion of the star on the MIPS detector. A reanalysis of our earlier 
observations of this system is consistent with our new result. 

Subject headings: infrared: planetary systems — planets and satellites: individual (y And b) - 
planetary systems — techniques: photometric — stars: individual (v And) 



1. INTRODUCTION 

The first thermal characterizations of highly irradiated 
extrasolar planetary atmospheres were made by mea- 
suring the flux decrement that occurs during secondary 
eclipse, when an extrasolar planet passes behind its host 
star (Deming et al. 2005; Charbonneau et al. 2005). 
This decrement gives an estimate of the (hemisphere- 
averaged) temperature of a planet's star-facing side at 
the time of eclipse and provides insight into the energy 
budgets of these hot worlds. Secondary eclipse (or, oc- 
cultation) observations have been widely interpreted as 
indicating two types of planetary atmospheres, namely 
planets with and without high-altitude temperature in- 
versions (Burrows et al. 2008; Fortney et al. 2008). How- 
ever, the number of free parameters in current models is 
such that in many cases it is still difficult to place strong 
constraints on a planet's atmospheric structure with the 
few data points available for most systems (Madhusud- 
han & Seager 2009). 
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Burrows et al. (2008) and Fortney et al. (2008) sug- 
gested that sufficiently high levels of irradiation pre- 
vent strong optically absorbing species from condens- 
ing and "raining out" of the upper atmosphere of hot 
Jupiters, thus directly linking high levels of incident 
stellar flux to the presence of a temperature inversion. 
However, subsequent secondary eclipse measurements 
with Spitzer/IRAC have complicated the picture and a 
straightforward connection between irradiation and in- 
versions now seems untenable. For example, TrES-3b re- 
ceives substantially more flux that does HD 209458b, yet 
the latter has an atmospheric inversion (Knutson et al. 
2008) while the former does not (Fressin et al. 2010). 
Thus, planetary classification will require more subtlety 
than a simple critical-flux level model can provide. Knut- 
son et al. (2010) have recently suggested a correlation 
between stellar activity and the absence of a tempera- 
ture inversion: in this hypothesis high-altitude absorbing 
species are photodissociated by the ultraviolet flux from 
an active star. It remains to be seen how this theory 
addresses the issue of temporally variable stellar activity 
(Shkolnik et al. 2008). 

Phase curves provide complementary insights into 
planetary atmospheres. If a system's total (star plus 
planet) infrared flux varies periodically and in phase 
with the planet's orbit, the variation can be attributed 
to spatially nonuniform radiation emitted by the planet. 
Such measurements have the potential to constrain the 
planet's circulation and heat redistribution patterns. If 
incident stellar flux were instantaneously re-radiated by 
the planet, the hottest region on the planet would be 
at the substellar point; such a phase curve is said to 
have zero phase offset. Nonzero phase offsets thus imply 
heat transport around the planet; for example, by advec- 
tion of absorbed stellar energy by global winds (Show- 
man et al. 2009), or by heating induced by atmospheric 
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gravity waves (Watkins & Cho 2010). Interpreting phase 
curves can be challenging because the brightest atmo- 
spheric region will also depend on the opacity structure 
of the atmosphere and the wavelength at which one ob- 
serves the system. 

The extrasolar planet v Andromedae b (v And b) 
is the first exoplanet for which a phase curve was re- 
ported. Harrington et al. (2006; hereafter H06) used 
Spitzer/MIPS to measure the 24 /im system flux at five 
epochs over one orbit and reported a finite amplitude 
phase curve consistent with zero phase offset, though 
here we report a new analysis and interpretation of our 
H06 data based on a better understanding of MIPS 
systematics. The 8 /tm observations of Cowan et al. 

(2007) also found variations with zero phase offset for 
HD 179949b. Observations of the less intensely irradi- 
ated planet HD 189733b at 8 /mi and 24 /im (Knutson 
et al. 2007, 2009) revealed a relatively small tempera- 
ture contrast between the planet's day and night sides 
and a 30° — 40° phase offset, indicating a moderate level 
of eastward heat redistribution from the warm dayside 
to the cool night side. It is important to note a pos- 
sible observational bias: the first two phase curves were 
sparsely sampled and are for non-transiting systems with 
unknown orbital inclinations. Though simulations sug- 
gest inclination should not substantially affect a planet's 
observed phase offset, the flux amplitude will be directly 
affected by inclination (Rauschcr et al. 2008); further- 
more, phase curve interpretations are more ambiguous 
without the absolute calibration provided by a secondary 
eclipse (Burrows et al. 2008). 

Several groups have hypothesized a connection be- 
tween temperature inversions and the magnitude of a 
system's phase offset (e.g., Burrows et al. 2008; Fortney 
et al. 2008). The favored (though unproven) cause of 
inversions is a species residing at high altitude that ab- 
sorbs optically but is transparent to infrared radiation. 
In this scenario an inverted atmosphere absorbs stellar 
energy at lower pressures where it should quickly rera- 
diate to space; in a non-inverted atmosphere the energy 
is absorbed much deeper, where it may circulate farther 
around the planet and cause a measurable phase offset. 
Showman et al. (2009) do a fair job of reproducing the 
HD 189733b phase curves, but they predict secondary 
eclipse depths for the more highly irradiated HD 209458b 
that do not match the observations of Knutson et al. 

(2008) ; Burrows et al. (2010) also model HD 209458b 
and they, too, do not match the observed eclipse depths 
especially well. Thus, our current understanding of the 
atmospheric structure and dynamics of even the best- 
characterized planets still appears to be incomplete. 

This is the context in which we obtained the high- 
cadence 24 /im phase curve of v And b described below. 
We introduce the v Andromedae system, and discuss our 
observations and data analysis, in Sec. 2. In Sec. 3 we 
describe the planetary temperature contrast and heat re- 
distribution implied by our analysis. We discuss possible 
interpretations of our results in Sec. 4, and conclude in 
Sec. 5. 



2. OBSERVATIONS AND ANALYSIS 
2.1. The v Andromedae System 



The planet v And b was one of the earliest reported 
hot Jupiters (Butler et al. 1997), and the three-planet 
v And system has been observed numerous times in the 
years since (Butler et al. 2006; Naef et al. 2004; Witten- 
myer et al. 2007; McArthur et al. 2010). The host star 
has a spectroscopic effective temperature of 6212 ± 64 K 
(Santos et al. 2004) and a directly-measured diameter of 
1.631 ±0.014 R Q (Baines et al. 2008). Spectroscopic and 
isochronal mass estimates generally agree on a mass of 
-1.3 M Q (Fuhrmann et al. 1998; Ford et al. 1999; Valenti 
& Fischer 2005; Takeda et al. 2007). Because this sys- 
tem's planets do not transit we do not know their physi- 
cal sizes; however, if we assume v And b is a typical hot 
Jupiter we can estimate its radius to be — 1.3 Rj s . 

Using a combination of radial velocity and astrome- 
try McArthur et al. (2010) recently determined the or- 
bits of the second and third planets in the v And sys- 
tem to be mutually inclined by 30°. They suggest this 
nonplanar system may result from planet-planet scatter- 
ing that could also have moved the innermost planet, 
v And b, into its current orbit at 0.059 AU. The small 
stellar reflex motion induced by v And b precluded a 
direct measurement of its orbital inclination, but their 
preliminary numerical simulations extending 10 5 yr sug- 
gest v And b's inclination may lie in the range —20-45° 
(implying a planetary mass of 2-3 Mj). Though their 
simulations did not fully explore the available parame- 
ter space, the inclination range McArthur et al. (2010) 
suggest for v And b is broadly consistent with the con- 
straints we place on its inclination in Sec. 3. 

2.2. Observations 

We observed the v And system with Spitzer's MIPS 
24 /im channel (Rieke et al. 2004) with observations 
spread across 1.2 orbits of v And b (—5 days) during 
February 2009. The observations consist of seven brief 
(—3000 seconds on target) observational epochs and one 
long, near-continuous observation —28 hours in length 
and centered at phase 0.5 (secondary eclipse for tran- 
siting systems in circular orbits, and the predicted time 
of flux maximum based on H06). Our integrations total 
18.5 hours. Spitzer breaks up observations into blocks 
of time called astronomical observation requests (AORs) 
for instrument scheduling purposes: our short observa- 
tions consist of three sequential AORs, and the long ob- 
serving sequence consists of 71 AORs. Altogether our 
data consist of 25 488 frames, each with an integration 
time of 1.57 seconds. We also observe a flux calibrator 
star, HD 9712, in three two-AOR epochs, for a total of 
1.3 hours of integration. The observations of v And by 
H06, which we reanalyze, consist of five AORs spaced 
over —5 days. 

2.3. Data Reduction 

We use the basic calibrated data (BCD) files gen- 
erated by version 18.14 of the MIPS data reduction 
pipeline (Masci et al. 2005). During MIPS observations 
the instrument and spacecraft dither the target star be- 
tween fourteen positions on the detector (SSC 2007, Sec- 
tion 8.2.1.2). As noted previously (Deming et al. 2005; 
H06; Knutson et al. 2009) the MIPS detector response 

8 Taken from the Extrasolar Planets Encyclopedia at http : / / 
exoplanet . eu 
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is spatially nonuniform and so we treat the observations 
as consisting of fourteen separate time series, modeling 
their systematics separately in the final fit. In each frame 
we measure the system flux and the position of the star 
on the detector by fitting a 100 x supersampled model 
MIPS PSF 9 generated using a 6200 K blackbody spec- 
trum. We shift and scale the model PSF to determine 
the best-fitting combination of background, stellar flux, 
and PSF position. Using position-dependent model PSFs 
does not significantly change our results, so in all our 
photometry we use a single model PSF generated at the 
center of the MIPS field of view. We exclude hot pixels 
from the PSF fit by setting to zero the weight of any 
frame's pixel that is more than 5c discrepant from the 
median value of that pixel for all frames taken at that 
particular dither position, and also exclude bad pixels 
flagged by the MIPS reduction pipeline. 

In each set of ^40 frames, the first frame has ^3% 
higher stellar flux, and the first several frames have a 
lower background, than the rest of the frames in the set. 
These effects may be related to the MIPS "first-frame" 
chip reset effect, though the effect we see is qualitatively 
different from the one described in the MIPS Instrument 
Handbook (SSC 2010). We exclude the first frame in 
each set from the remainder of our analysis, removing 
708 frames. We further exclude the first three contiguous 
AORs (700 frames, ~ 40 minutes), which are markedly 
discrepant from the final time series. These first data 
were taken soon after a thermal anneal of the MIPS de- 
tector, and during these observations we see anomalous 
readings in the 24^tm detector anneal current, the scan 
mirror temperature, and the MIPS B side temperature 
sensor. 

The initial extracted photometry reveals a clear sinu- 
soidal flux variation, but with additional variability cor- 
related with the subpixel motion of the star on the de- 
tector, as shown in Figures 1 and 2. This effect is likely 
from imperfect flat fielding rather than an intrapixel sen- 
sitivity variation as seen in the Infrared Array Camera's 
3.6 and 4.5 /Ltm channels (Reach et al. 2005; Charbonncau 
ct al. 2005), because several time series with similar in- 
trapixel locations exhibit anticorrelated flux variations. 
We tried creating a flat field by median-stacking all the 
BCD frames after masking the region containing the tar- 
get star, but applying this flat field to the data before 
computing photometry did not reduce the amplitude of 
the position-correlated photometric variability. As we 
describe below, we suspect intermittent, low levels of 
scattered light interfere with our ability to construct a 
sufficiently accurate flat field 

Instead, we treat this variability by removing a lin- 
ear function of position from the computed photome- 
try at each of the fourteen dither positions. Fitting 
functions with a higher-order dependence on position, 
or including cross terms, does not change our final re- 
sults and increases the Bayesian Information Criterion 
(BIC = x 2 + k In N, where k is the number of free pa- 
rameters and N is the number of points). We remove the 
systematic effects, and assess possible phase functions, by 

fitting the j^ 1 data point at the i^ 1 dither position (fij) 

9 Generated using Tiny Tim; available at http: //ssc . spitzer . 
caltech . edu/ 



with the following equation: 
fij = (a + frsin [^ orD i - $o]) (1 + Cj + ckxij + e^yy) 

This equation contains three astrophysical parameters 
and 42 instrumental parameters that account for the 
nonuniform detector response. The astrophysical pa- 
rameters of interest are the average system flux a and a 
time (i)-dependent sinusoidal phase function with known 
planetary orbital frequency ^ orD but unknown ampli- 
tude b and phase offset <f>o- The remaining parameters 
represent systematic effects to be removed: residual sen- 
sitivity corrections Ci and linear dependence on detector 
position (di,ei) for each of the fourteen time series. 

To prevent parameters a and b from floating we artifi- 
cially set ci to satisfy the relation Hj(l + q) = 1. Ignor- 
ing the (di,ei) factors does not change our final result 
for the 2009 dataset but increases the BIC, indicating a 
poorer fit to the data. For the purposes of backwards 
comparison, we also apply our analysis to the observa- 
tions of H06. Due to the limited temporal coverage of 
this dataset including the (dj,ej) does not improve the 
fit; therefore we use only 17 parameters when reanalyz- 
ing this earlier dataset. We otherwise apply the same 
data reduction steps as described above. 

In this work we consider only a sinusoidal phase func- 
tion. Although we recognize that the phase curve will 
not be truly sinusoidal in shape, such a model is simple 
to work with and has a straightforward interpretation 
as a two-hemisphere "orange slice" model (Cowan et al. 
2007; Cowan & Agol 2008, see also Sec. 3.2 below). Since 
the motion-correlated flux variation only depends on rel- 
ative motion (rather than on absolute detector position) 
we normalized the (x^-, y^ ) in Eq. (1) by subtracting the 
mean position in each of the fourteen time series. 

2.4. Calibration and instrument stability 

Our continuous photometry reveals that the MIPS 
background flux changes discontinuously from one AOR 
to the next. Stellar photometry is not similarly affected. 
The background flux varies at the level of 0.1%, com- 
parable in amplitude to the expected planetary signal. 
We also see these discontinuities in background flux dur- 
ing MIPS observations of HD 189733b (Knutson et al. 
2009) and HD 209458b (unpublished; Spitzer Program 
ID 40280). 

Because the measured background level varies with the 
Spitzer AORs it is extremely unlikely that this variability 
is of astrophysical origin. It is also unlikely that the back- 
ground variability results from the calibration process 
because we see the same effect in both the raw and cali- 
brated data products. We observe no correlation between 
the background variability and the various reported in- 
strumental parameters, though we cannot rule out either 
intermittent scatter from other sources or slight changes 
in the detector bias. A global sensitivity drift does not 
seem to be the culprit because the changes in background 
flux are uncorrelated with stellar photometry. 

Although we were unable to determine the cause of 
the background variations, we suspect that they are due 
to small changes in the amount of scattered light reach- 
ing the MIPS detector. Variable scattered light could 
also explain our inability to remove the motion-correlated 
photometric jitter with an empirical flat field. Stel- 
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lar photometry was substantially more stable on short 
timescales than was the background in all of the extended 
MIPS 24/im observations we examined, so we use this 
photometry in our subsequent analysis. 

Our observations include a flux calibrator star to check 
MIPS's long-term photometric stability. MIPS 24 ^m 
photometry is known to be stable at the 0.4% level over 
several years (Engelbracht et al. 2007); if it were this 
variable on short timescales we would be unable to dis- 
cern the expected planetary emission. We observed the 
Kl III star HD 9712, taken from a catalog of bright inter- 
ferometric calibrator stars (Merand et al. 2005), during 
six AORs. In our reduction of these data we do not apply 
a correction for pixel motions and we achieve a repeata- 
bility of ~ 10~ 4 . These observations, shown as the red 
triangles in Figure 3, imply that the MIPS 24 fim detec- 
tor was stable over the course of our observations, and 
so we rule out detector sensitivity drifts as the source of 
the observed flux variations. 

The level of MIPS photometric stability was not well 
known early in the Spitzer mission, so H06 used the back- 
ground, attributed to smoothly-varying zodiacal emis- 
sion, as a calibrator to adjust the photometry (see H06 
Figure IB). As we now suspect the MIPS background 
variations to result from scattered light in the instru- 
ment, there is no longer any need, nor possibility, to use 
the zodiacal light as a calibrator. The H06 data were 
part of a preliminary Spitzer program to assess variabil- 
ity of several systems for subsequent study in programs 
such as ours. Its five v And AORs were taken at separate 
epochs and there was no calibrator star, so it was impos- 
sible to make the assessment described above. We have 
reanalyzed the H06 data using our new procedure and 
find values consistent with those plotted in H06 Figure 
1A (see our Figure 3) and also with our new data. Thus, 
H06 did still present the first orbital phase variation for 
an exoplanet. The phase curve fit here to the H06 data 
provides accurate parameters for that dataset. 

3. RESULTS 

3.1. System flux 

Although our primary science result - the planetary 
phase curve described below - is inherently a relative 
measurement, our observations also allow us to measure 
precise, absolute 24/mi photometry for the v And sys- 
tem. We measure F„ = 0.488 Jy and 0.490 Jy for the 
2009 and 2006 datasets, respectively. These fluxes dif- 
fer by 0.4%, which is at the limit of the MIPS 24^m 
precision; we therefore report the mean system flux as 
0.489 ± 0.002 Jy. This value is significantly discrepant 
from the IRAS 25/zm flux of 0.73±0.05 Jy (Moshir 1989), 
but it is consistent with a Kurucz (1979) stellar spectrum 
tied to optical and near-infrared photometry of v And 
from the Tycho-2 (H0g et al. 2000) and 2MASS (Skrut- 
skie et al. 2006) catalogs. 

3.2. A two -hemisphere model 

As shown in Figure 3 and discussed below, we de- 
tect a flux variation consistent with the orbital period 
of v And b. This measurement allows us to put tight 
constraints on the temperature contrast and phase offset 
of the planet. We interpret the observed flux variation 
in the context of a planet composed of two blackbody 



hemispheres of constant temperature - i.e., a two- wedge 
"orange slice" model (c.f. Cowan et al. 2007). Suffi- 
ciently precise observations of such a bifurcated planet 
will reveal a flux variation with peak-to-trough ampli- 
tude 
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which normalizes the full amplitude of the planetary flux 
variation AFp by the mean flux from the system (F) ; in 
Eq. (1) AFp /(F) = 2b/ a. The quantity 7 accounts for 
the star being fainter in the mid-infrared than a black- 
body with temperature T^ff! at 24/zm we set 7 = 0.8 
based on the models of Kurucz (1979). Thus, measuring 
AFp/ (F) gives the hemispheres' brightness temperature 
contrast relative to the stellar flux, modulo a sin i ambi- 
guity. By assuming a planetary albedo Ab and assuming 
all emitted radiation is reprocessed starlight, a second 
(bolometric) relation obtains: 



-pi- 



(3) 



Thus, we implicitly assume that each hemisphere of 
our model planet emits a bolometric flux equal to that 
of a blackbody with the 24 /im brightness temperature in 
Eq. (2). Subject to this assumption and given known or 
assumed values for the albedo Ab and planetary radius 
Rp, one can use Eqs. (2) and (3) to determine the hemi- 
spheric temperatures Tp± and T P i at arbitrary orbital 
inclinations. 

Setting Tp2 = and T P \ = Tpi m &x gives the min- 
imum planetary radius capable of reproducing the ob- 
served flux variation, AFp / (F) , as a function of the or- 
bital inclination. Because this relation depends on incli- 
nation as (sini) , measuring AFp /(F) gives an upper 
limit to the planet's surface gravity: 



< 



G(m sin i) B v (Tpi >m ax) / AF P 
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(4) 



3.3. Model fits 

We determine the best-fit parameters using Powell's 
(1964) method for multivariate minimization (the SciPy 
function optimize. fmin_powell), and assess their uncer- 
tainties and correlations with a Metropolis-Hastings, 
Markov-chain Monte Carlo analysis (MCMC; see Press 
et al. 2007, , section 15.8). Table 1 (for the 2009 data) 
and Table 2 (for the 2006 data) report these results. Ta- 
ble 3 lists the astrophysical parameters of interest. We 
list the x 2 an d BIC values for the fits in Table 4. Pa- 
rameter uncertainties are estimated from distributions 
generated using the kernel density method (KDE, im- 
plemented using the SciPy function stats. gaussian_kde) 
by determining the parameter values with equal KDE 
frequency that enclose 68% of the distribution. 

MCMC analysis evolves an initial set of parameters in 
a way that is ultimately representative of their underly- 
ing probability distributions. For our Markov chain we 
choose a step size to give approximately a 30% step ac- 
ceptance rate. To adequately sample the full parameter 
space we found it necessary to run the Markov chains 
longer for the 2006 dataset than for the 2009 dataset. 
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For the 2009 dataset we first ran the chain for f0 6 burn- 
in steps and discarded these; we then ran the chain for 

2 x 10 7 steps, saving every 1000^ step. For the 2006 
dataset our procedure is the same but the burn-in phase 
lasted for 10 7 steps and the chain was then run for 5 x 10 7 

steps, saving every 1000^ step. We inspected correlation 
plots for all possible parameter pairs in both analyses to 
ensure adequate coverage of phase space and to assess 
parameter correlations. 

All the one-dimensional parameter distributions are 
unimodal and approximately Gaussian in shape. We see 
some correlations between the di and , which is unsur- 
prising given the degree of correlation between the X and 
Y components of motion as shown in the middle panels 
of Figure 1. More surprising is a correlation between the 
mean system flux, a, and the phase offset, $o, as shown in 
Figure 4. Using a simulated dataset with white noise, we 
confirmed that when forcing a fit to a sinusoid of known 
period, a slightly higher mean value be counteracted by 
a slightly lower phase offset; however, we observe the op- 
posite correlation. In any case no significant correlation 
is apparent between the phase curve amplitude AFp/(F) 
and the phase offset, which arc the primary quantities of 
interest for our analysis. The best \ 2 from the MCMC is 
consistent to within a small fraction of the uncertainties 
with the optimizer values. 

The model in Eq. (1) provides a good fit to both the 
astrophysical flux modulation and the instrumental flux 
variations at each of the fourteen dither positions, as 
shown in Figure 2. We plot the photometry after removal 
of the systematic effects in Figure 3, along with the best- 
fitting sinusoidal phase curve, for both the 2006 and 2009 
datasets. Both datasets appear to vary approximately in 
phase; this coherence is a strong argument that the flux 
modulation we see is due to the planet v And b. The 
goodncss-of-fit statistics (% 2 and BIC) for both datasets 
are listed in Table 4. The high \ 2 for the 2009 dataset 
probably results from a somewhat non-sinusoidal phase 
curve and from the residual systematics apparent in Fig- 
ure 3. 

3.4. Phase curve amplitude 

We measure values of AF P /(F) of AF P /(F) = 
0.001300 ± 0.000074 for the 2009 data and 0.00090 ± 
0.00022 for our reanalysis of the 2006 data, as shown in 
Table 3. Using the absolute calibration from Sec. 3.1, 
we find absolute peak-to-trough phase curve amplitudes 
of 0.636 ± 0.036 mJy and 0.44 ± 0.11 mJy for the 2009 
and 2006 datasets, respectively. Thus the detection of 
the phase curve amplitude at both epochs is statistically 
significant at the > 4cr level and is substantially smaller 
than reported by H06 (due to the calibration issues dis- 
cussed above). The two epochs' phase curve amplitudes 
are consistent at the 1.7a level; thus there is no evidence 
that the planetary emission exhibits inter-epoch variabil- 
ity. The lack of variability is consistent with the recent 
results of Agol et al. (2010), who set an upper limit of 
2.7% on HD 189733b's dayside flux variations. 

As v And b's orbit is inclined toward face-on, a greater 
intrinsic temperature contrast is required to generate the 
observed flux variation. Using Eq. (2) and our measure- 
ment of the phase curve amplitude we determine the ex- 
pected day/night contrast ratio and plot it in the up- 



per panel of Figure 5. We also plot the upper limits 
on the day/night contrast assuming planetary radii of 
1.3 (1.8) Rj; the implication is that the planet's orbital 
inclination angle is likely > 28° (14°). These limits com- 
plement the preliminary limits on v And b's orbital in- 
clination from the stability modeling of McArthur et al. 
(2010), which suggest i ~ 20° - 45°. 

Invoking Eq. (3) allows us to determine the brightness 
temperatures of the planetary hemispheres in our model 
at each inclination angle. We assume zero albedo (c.f. 
Rowe et al. 2008) and a planetary radius of 1.3 Rj and 
plot the hemispheres' temperatures and 3cr limits in the 
lower panel of Figure 5. This sets a lower bound to the 
temperature contrast between the two hemispheres to be 
Tpi — Tp2 > 900 K. The hotter hemisphere's tempera- 
ture remains in the range ^1700-1900 K as we vary the 
radius from 1.0 Rj to 1.8 Rj] however larger radii re- 
sult in higher temperatures for the cooler hemisphere to 
maintain the measured flux and thus decrease the tem- 
perature contrast. 

Using Eq. (4) we find that v And b's surface gravity 
is < 2100 cm s -2 with 3er confidence: this result is in- 
dependent of assumptions about the planet's radius or 
orbital inclination. For a hot hemisphere temperature of 
^1800 K (c.f. Figure 5) this limit on the surface grav- 
ity implies an atmospheric scale height > 300 km. In 
Figure 6 we plot the allowed regions of mass-radius pa- 
rameter space against the known population of transit- 
ing extrasolar planets. Thus, our measurements suggest 
that v And b has a lower surface gravity than Jupiter, 
HD 189733b, and a number of other transiting extraso- 
lar planets. This result demonstrates that v And b is 
indeed a gaseous Jovian planet, but we cannot deter- 
mine whether it is a highly inflated planet or whether it 
is dominated by a sizeable rocky core. 

3.5. Phase offset 

Because v And b does not transit its host star we know 
its orbital ephemeris less precisely than we do for transit- 
ing planets. We used the Systemic Console 10 (Mcschiari 
et al. 2009) to reanalyze the published radial velocity 
data of v And (Butler et al. 1997; Naef et al. 2004; But- 
ler et al. 2006; Wittenmyer et al. 2007; McArthur et al. 
2010) using Systemic's Levenberg-Marquardt algorithm 
and ignoring system stability constraints. We obtain or- 
bital parameters consistent with those of McArthur et al. 
(2010). By providing the covariances between the vari- 
ous fit parameters, this reanalysis allows a substantially 
more precise estimate of the planetary ephemeris than is 
available from the literature. We compute a time of zero 
relative radial velocity (phase 0.5, or secondary eclipse in 
a circular transiting system) of JD = 2 454 868.78 ± 0.07 
(ler), which corresponds to an uncertainty of ~ 6° in 
determining the phase offset. 

Assuming a two-hemisphere model, we find a phase off- 
set of 84.5° ± 2.3° relative to the computed ephemeris. 
This uncertainty may be an underestimate since we arti- 
ficially constrain our phase curve to be sinusoidal. Knut- 
son et al. (2009) discuss the artificially low uncertain- 
ties obtained from fitting to an arbitrarily chosen model, 
though here we have substantially broader phase cover- 
age than was available in that study. The system flux 

10 available from http://oklo.org/downloadable-console/ 
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reaches a maximum before phase 0.5, indicating that the 
brighter hemisphere is offset to the east of the substellar 
point, as observed for HD 189733b (Knutson et al. 2007, 
2009; Agol et al. 2010). 

This large phase offset is strikingly different from the 
near-zero phase offset reported by H06; this difference 
is due to the choice of system calibration as discussed 
above. From our reanalysis of the 2006 data we find 
a phase offset of 57° ± 21° relative to our ephemeris. 
The phase offsets at the two epochs are consistent at the 
1.3cr level, and thus there is no evidence for inter-epoch 
variability in the phase offset. 

4. DISCUSSION 

The most striking result of our analysis is the large 
phase offset evident in the light curve. The direc- 
tion of the phase offset is broadly consistent with the 
prediction by many circulation models of a large-scale, 
high-velocity, eastward-flowing jet on hot Jupiters (Cho 
et al. 2003; Cooper & Showman 2005; Cho et al. 2008; 
Showman et al. 2009; Burrows et al. 2010; Rauscher & 
Menou 2010; Thrastarson & Cho 2010) and as seen on 
HD 189733b (Knutson et al. 2007, 2009). However, the 
magnitude of the phase offset is far larger than is pre- 
dicted at the low pressures characteristic of the 24 fim 
photosphere (Showman et al. 2009). 

A partial explanation for such a large phase offset could 
be that v And b's atmosphere is substantially transpar- 
ent to the incident stellar flux, with the result that the 
insolation is deposited at sufficient depth for substantial 
advection to occur. In this case, the planet might lack 
an atmospheric temperature inversion. Indeed, the large 
phase offset is more consistent with expectations for a 
planet lacking a temperature inversion than with expec- 
tations for a planet with an inverted atmosphere (c.f. 
HD 189733b vs. HD 209458b; Showman et al. 2009). 
Cooper & Showman (2005, 2006) predicted maximum 
hemisphere-averaged temperatures to be offset by > 60° 
for an irradiated planet with no high- altitude absorbers, 
but it is unlikely that this would translate into high phase 
offsets at 24 fim due to the high altitude of the photo- 
sphere at this wavelength. The absence of a temperature 
inversion is at odds with the recent proposal by Knutson 
et al. (2010) that, due to the lack of activity apparent 
in optical spectra of v Andonc would expect v And b to 
have a temperature inversion. Thrastarson & Cho (2010) 
also show clear shifted hot regions at higher pressure; the 
deep vortices produced by their simulations are several 
hundreds of degrees hotter than the surroundings. 

Setting aside the phase offset for a moment, the phase 
curve amplitude we measure could be directly interpreted 
in the context of the models of Showman et al. (2009) as 
a typical phase curve for a planet with orbital inclina- 
tion of 50°-60°. However, simulations with greater phase 
offsets tend to show lower phase curve amplitudes (e.g. 
Showman et al. 2009; Burrows et al. 2010); thus it is 
difficult to simultaneously explain both amplitude and 
offset and it is difficult to reconcile current theory with 
our observations. 

Alternatively, the large phase offset we see could rep- 
resent reradiation of thermal energy deposited in shock 
fronts in the planet's atmosphere. Many simulations 
predict supersonic equatorial jets on hot Jupiters that 
carry substantial kinetic energy; shocks could manifest 



themselves where the jet transitions to subsonic speeds. 
Rauscher & Menou (2010) note that their simulations, as 
well as those of Showman et al. (2009), exhibit structures 
which could be interpreted as shocks - however, these 
models do not explicitly treat shock physics. Dobbs- 
Dixon et al. (2010) observe similar features in their mod- 
els using an additional artificial viscosity factor to simu- 
late shock behavior, and note that in high-altitude, high- 
velocity regions the energy carried by kinetic energy be- 
comes comparable to the enthalpic energy. Watkins & 
Cho (2010) have recently suggested that gravity waves 
in the atmosphere of a hot Jupiter can also heat the 
planet's upper atmosphere. It is unclear whether either 
shocks or gravity waves can deposit sufficient energy in 
the 24 fim photosphere to cause the large phase offset we 
see. We look forward to further research into these top- 
ics to determine whether these or other phenomena can 
explain our observations. 

We must also consider the possibility that the peri- 
odic flux modulation we see is intrinsic to the star rather 
than emanating from the planet. Shkolnik et al. (2005, 
2008) report evidence for intermittently periodic stellar 
activity in this system correlated with v And b's or- 
bital period. They interpret this periodicity as a possi- 
ble magnetospheric star-planet interaction, but they de- 
tect this periodicity at only some of their observational 
epochs. The consistency between our analysis of the 
2006 and 2009 MIPS data suggests we are not seeing such 
a transient phenomenon, although analysis of additional 
activity measurements (e.g., from the Keck Observatory 
Archive) spanning the time of our Spitzer observations 
would help solidify this claim. Shkolnik et al. (2005) also 
report variations consistent with zero phase offset, which 
disagree with our observed phase offset. In addition, 
24 fim stellar variability at the level we observe would 
imply much greater variability at optical wavelengths: 
this is not observed (Henry et al. 2000). Finally, as dis- 
cussed by H06, energy considerations indicate that such 
an intense star-planet interaction would cause the planet 
to spiral into the star in < 10 7 yr. Thus, the 24 fim flux 
variations are likely to be of planetary origin. 

Although a general framework exists within which to 
interpret observations of exoplanetary atmospheres, our 
understanding is still extremely limited. Madhusudhan 
<fc Seager (2009) demonstrate that some planets with 
claimed inversions can also be fit by non-inverted at- 
mospheric models due to much greater number of free 
model parameters versus the limited number of obser- 
vational constraints. In addition, there may be rea- 
son to question the reliability of some of the circulation 
models currently in use. Thrastarson & Cho (2010), in 
their extensive exploration of initial condition parame- 
ter space, have recently shown the extreme susceptibility 
of at least one circulation model to minute variations 
in initial conditions, resulting in substantial variability 
in final "steady-state" temperature contrasts and phase 
offsets. This chaotic behavior calls into question the 
ability of at least some models to make accurate, qual- 
itative predictions about any of the quantities we are 
interested in. Other models (e.g., Dobbs-Dixon et al. 
2010) exclude planets' polar regions from their simula- 
tions, where large vortices are often seen to form (Cho 
et al. 2003, 2008; Rauscher et al. 2008). While the polar 
regions present a small cross-sectional area in transiting 
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systems, for non-transiting systems such as v And b po- 
lar emission will constitute a larger component of the ob- 
served system flux. A fully consistent three-dimensional 
circulation geometry is essential for comparison to our 
observations. Nevertheless, it is important to remember 
that general circulation models have trouble predicting 
global weather patterns even for relatively well-studied 
solar system planets; thus, a comprehensive, quantitative 
understanding of extrasolar planetary atmospheric dy- 
namics will likely remain elusive for some time to come. 

5. CONCLUSIONS AND FUTURE WORK 

We have described a new 24 /im phase curve, which we 
interpret as being due to emission from the planet v An- 
dromedae b modulated by the planet's orbit. Using a 
simple two-hemisphere model we determine the peak-to- 
trough phase curve amplitude to be 0.001300 ±0.000074. 
This result suggests an average "hot side" temperature 
of ~1800 K; for an average-sized hot Jupiter (1.3 Rj) 
this implies a hemisphere-averaged planetary temper- 
ature contrast of ^900 K and an orbital inclination 
i > 28°. 

We find a phase offset of 84° ± 2° ± 6° , where we break 
the uncertainty into the error relative to our ephemeris 
and the error in our ephemeris, respectively. Such a large 
phase offset is difficult to reconcile with most current 
models. The phase curve is hotter hemisphere is off- 
set to the east, as previously observed for HD 189733b 
(Knutson ct al. 2007, 2009). 

We reanalyze our earlier (H06) observations of this 
system and find a phase curve amplitude of 0.00090 ± 
0.00022 and a phase offset of 57° ± 21° ± 6°. There is 
no evidence for inter-epoch variability in the planetary 
phase curve. This is primarily due to the large uncer- 
tainties from the 2006 dataset and demonstrates the diffi- 
culty in measuring such variability with sparsely-sampled 
phase curves. 

There are substantial challenges in interpreting a phase 
curve observed at only a single wavelength due to the 
degeneracies between planetary radius, orbital inclina- 
tion, and atmospheric composition and structure (Bur- 
rows et al. 2008). Some of these difficulties could be 
mitigated with phase curve measurements at additional 
wavelengths - ideally from space (i.e., warm Spitzer) but 
also potentially from ground-based near-infrared obser- 
vations (e.g., Barnes ct al. 2010). For example, com- 
mensurate phase curve amplitudes at both 3.6 /im and 
4.5 /im would suggest that v And b does indeed lack 
a temperature inversion. On the other hand, differing 
phase curve amplitudes at these wavelengths could sug- 
gest an inversion and be more difficult to reconcile with 
the large phase offset at 24 /im. Phase offsets in these or 
other scenarios would also depend on the particular at- 
mospheric temperature structure of the planet. Phase 
curves at 24 /im currently exist for only two planets, 
v And b and HD 189733b (though unpublished data ex- 
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ist for HD 209458b). Cool, inversionless HD 189733b 
is probably the best-characterized extrasolar planet, and 
additional phase curves are already being observed for 
this object at 3.6 /im and 4.5 /im. Phase curves of ad- 
ditional planets at multiple wavelengths are essential to 
ensure that our evolving views of the atmospheres of hot 
Jupiters are not biased by unbalanced datasets. 

Whatever the cause of the substantial energy trans- 
port implied by the large phase offset we measure, if 
this phenomenon occurs in other (transiting) systems 
there are important implications for transmission spec- 
troscopy. Line-of-sight effects cause optical transmission 
spectra to probe pressures comparable to those probed 
by mid-infrared emission; temperatures near the plane- 
tary terminator of ^1800 K (as we observe) should easily 
be detectable with ground- or space-based spectra. Sing 
et al. (2008) and Lecavelier des Etangs et al. (2008) infer 
a terminator temperature on HD 209458b of 2200±260 K 
at 33 ± 5 mbar - roughly at the expected 24 /im photo- 
sphere - which is not dissimilar from the hot terminator- 
centered hemisphere we observe on v And b. Further 
phase curve and transmision spectra of additional sys- 
tems are needed to determine whether this hot, high- 
altitude terminator measurement results from a mecha- 
nism similar to what we observe on v And b. 

We note that v And b is too bright to be observed 
photometrically at shorter wavelengths with the planned 
James Webb Space Telescope (though spectroscopy may 
be feasible). Given JWST's magnitude limits and the 
expected high demand for its observing time, the com- 
munity should consider a dedicated space-based mid- 
infrared photometry and spectroscopy mission (Vasisht 
et al. 2008). Such a mission would allow uninterrupted 
long-term monitoring of nearby hot Jupiter systems. 
This would provide high-precision measurements of these 
systems' thermal emission and energy distributions, and 
possibly provide the first definitive evidence of dynam- 
ical meteorological processes - weather - on extrasolar 
planets. 
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the California Institute of Technology under a contract 
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residence at the Kavli Institute for Theoretical Physics, 
funded by the NSF through grant number PHY05-51164. 
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TABLE 1 
Phase curve and calibration 
parameters from eq. (1) (2009 data) a . 
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a Errors quoted are the 68.3% confidence 
limits. The first three parameters are 
the mean system flux a, the phase curve 
half-amplitude b, and the phase offset <I>o- 
The other coefficients represent fits to the 
nonuniform detector response at each of 
the fourteen MIPS dither positions. The 
d are dimensionless, and the di and ei are 
in units of pixel -1 . 

b The absolute accuracy of the mean sys- 
tem flux a is 0.4% as discussed by Engcl- 
bracht ct al. (2007). 

c <&o is measured relative to our computed 
ephemeris, as discussed in Sec. 3.5. 
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TABLE 2 

Phase curve and calibration 
parameters from eq. (1) (2006 

DATA) a . 
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a Errors quoted are the 68.3% confidence 
limits. The first three parameters are 
the mean system flux a, the phase curve 
half-amplitude b, and the phase offset c^o- 
The other coefficients represent fits to the 
nonuniform detector response at each of 
the fourteen MIPS dither positions. The 
Ci are dimensionless, and the di and ei are 
in units of pixel -1 . 

b The absolute accuracy of the mean sys- 
tem flux a is 0.4% as discussed by Engel- 
bracht et al. (2007). 

c $o i s measured relative to our computed 
ephemeris, as discussed in Sec. 3.5. 



TABLE 3 

ASTROPHYSICAL PARAMETERS OF INTEREST 





2006 


2009 


AF/F 


0.00090 ± 0.00022 


0.001300 ± 0.000074 


$0 


57° ± 21° 


84.5° ± 2.3° 



TABLE 4 

GOODNESS-OF-FIT STATISTICS 





2006 


2009 


x 2 


814.1 


38422.7 


h 


17 


45 


N 


838 


23884 


BIC 


928.6 


38876.4 
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Fig. 1. — Removal of pointing-correlated photometric variability. The top panel shows photometry after bringing all fourteen individual 
time series to a common level. A large-scale sinusoidal flux variation is evident, but so is a shorter-scale "ripple" (near phase 0.5); this 
ripple is correlated with the motion of v And on the MIPS detector, plotted as X and Y in the middle two panels. The bottom plot shows 
the final, cleaned photometry after removing this correlation. For display purposes the full dataset has been averaged over each AOR in 
this figure. Uncertainties on the X and Y points are typically ~10~ 3 pixel, while the uncertainties of the photometric points (not plotted 
for clarity) is typically ~1.5 X 10~ 4 . 
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Orbital Phase Orbital Phase 

Fig. 2. — Raw photometry (dots) and the best-fit model (solid lines) at each of the fourteen dither positions. The measured flux varies by 
several percent from one position to another, as evidenced by the different scales in each panel. The downward trend evident in all panels 
is due to the decreasing flux from the system, shown more clearly in Figure 3. For display purposes the data have been averaged over each 
AOR and we plot only the continuous-observing segment of our observations. The precision of individual points (not plotted for clarity) is 
~5x 10" 4 . 
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Fig. 3. — Phase curve of the v And system, phased to the orbit o innermost planet (black circles). Position-dependent sensitivity 
effects have been removed and for display purposes we have averaged the data over each AOR. The best-fit sinusoid (black dotted line) 
exhibits a phase offset of ~ 80°, consistent with a planetary "hot spot" advected almost to the planet's day-night terminator. The light 
green circles are our reanalysis of the data from Harrington ct al. (2006), and the light green dashed line shows the best-fit sinusoid to these 
data; the phase coherence between the 2006 and 2009 datasets is consistent with flux modulated by the innermost planet's orbit. The red 
triangles show our flux calibrator observations, which are consistent with a constant detector sensitivity. 



b [mJy] Phi_0 [rad] 

0.20 0.25 0.30 0.35 0.40 1.301.351.401.451.501.551.60 

I 1 1 1 1 1 1 1 1 1 1 1 1 1 0.48850 




Fig. 4. — Two-dimensional confidence intervals on the astrophysical parameters of interest in Eq. (1): the relative phase curve amplitude 
AFp/F, the absolute mean system flux a, the phase curve half-amplitude b, and the phase offset <E>o. The solid lines are the contours that 
enclose 68.3%, 95.4%, and 99.7% of the parameter space from the 2009 dataset. The 'x' in each panel marks the best-fit parameter listed 
in Tables 1 and 3. 
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Fig. 5. — Day-night temperature contrast as a function of orbital inclination angle, assuming a planet with zero albedo, (a) Our 
measurement of the phase curve amplitude and Eq. (2) constrain the allowed day/night contrast to lie between the dot-dashed 3cr limits. 
The maximum allowable contrast for a planet of 1.3 (1.8) Rj is shown as the dashed (dotted) lines, suggesting a lower limit on the 
inclination angle of 25° (15°). (b) Considering Eq. (3) and assuming a radius of 1.3-R,/ allows us to determine the temperature of both 
hemispheres as a function of inclination (solid lines, with dot-dashed 3cr limits). This radius implies a planetary temperature contrast of 
> 900 K. The cooler hemisphere is more sensitive to changes in planetary radius, though both temperatures increase as radius increases. 



New v And b Phase Curve 



15 




0,6 0.5 1.0 1.5 2.0 2.5 3.0 

Mass/Mjup 

Fig. 6. — Mass/radius constraints for v And b from Eq. (4). The shaded area is the portion of mass-radius space excluded at the 3cr 
level. Points are known transiting extrasolar planets; measurement uncertainties have been omitted for clarity. Extrasolar planets that 
have been observed with MIPS are indicated by a '+' (HD 189733b) and a 'x' (HD 209458b), while Jupiter and Saturn are marked with a 
'J' and 'S,' respectively. The dot-dashed lines represent lines of constant surface gravity with g = 10 3 cm s" " 2 and 3.2 X 10 3 cm s . Our 
model and phase curve measurement constrain g < 2100 cm s — 2 at the 3a level. 



